The Magnetoviscous-thermal Instability 
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ABSTRACT 



Accretion flows onto underluminous black holes, such as Sagittarius A* at 
the center of our galaxy, are dilute (mildly collisional to highly collisionless), 
optically thin, and radiatively inefficient. Therefore, the accretion properties of 
such dilute flows are expected to be modified by their large viscosities and ther- 
mal conductivities. Second, turbulence within these systems needs to transport 
angular momentum as well as thermal energy generated through gravitational 
infall outwards in order to allow accretion to occur. This is in contrast to clas- 
sical accretion flows, in which the energy generated through accretion down a 
gravitational well is locally radiated. In this paper, using an incompressible fluid 
treatment of an ionized gas, we expand on previous research by considering the 
stability properties of a magnetized rotating plasma wherein the thermal con- 
ductivity and viscosity are not negligible and may be dynamically important. 
We find a class of MHD instabilities that can transport angular momentum and 
thermal energy outwards. They are plausible candidates to describe accretion in 
radiatively inefficient accretion flows (RIAFs). We finish by discussing the impli- 
cations for analytic models and numerical MHD simulations of mildly dilute or 
collisionless astrophysical plasmas, and immediate directions for further research. 



1. Introduction 



The study of astrophysical rotating plasmas in general, and accreting plasmas in 
particular, has undergone a renaissance within the past twenty years. |Balbus fc Hawley" 



(1991); Hawley & Balbus (1991) reapplied the magnetorotational instability (MRI, 



Velikhov (1959); Chandrasekhar (I960)) into an appropriate astrophysical context. They 
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demonstrated that in rotating astrophysical plasmas, even arbitrarily small magnetic forces 
can provide a channel by which free energy source of outwardly decreasing angular velocity, 
rather than angular momentum, can destabilize a Rayleigh-stable flow, generate magnetic 
fields with energies of order the thermal energy, and drive the right form of turbulence that 



allows for accretion to occur. Numerous numerical simulations (Hawley et al. 1996 Wardle 



1999 Sano & Stone 2002 De Villiers & Hawley 2003 Fromang et al. 2004) have borne out 



the fact that magnetic fields can play an essential role in driving efficient accretion for a wide 
class of astrophysical objects. The nonlinear MRI is currently the most plausible candidate 
to explain high mass-rate, radiatively efficient accretion onto massive and supermassive 
central galactic black holes (M ^ 1O 7 M ), which manifests itself as active galactic nuclei 
(AGNs) and quasars. 

Although AGNs and quasars are some of the most powerful and most energetic compact 
objects in the universe, at recent epochs in the history of the universe (at redshifts less than 
1) they are rare. Evidence from Very Large Baseline Interferometry (VLBI) measurements 
of water masers, from gas spectroscopy, and from stellar dynamics give strong evidence for 
the existence of massive and supermassive black holes in galaxies that possess a central 



bulge (see Richstone et al. (1998) for an exhaustive review on diagnostic tools for the 
detection of central galactic black holes). However, recent X-ray surveys of galactic and 



extragalactic sources (Tueller et al. 2008, 2010) has demonstrated that there is strong 



observational evidence of approximately 230 local AGNs, located within Seyfert galaxies, 
within 400 million light-years. The relative transparency of active galactic nuclei hard X-ray 



emission to the high column density material obscuring them (Treister & Urry 2005), and 



the lack of confusion with other sources of emission in the AGNs host galaxy (Markwardt 



et al. 2005), makes hard X-ray emission an optimal spectral regime for unbiased AGN 
searches. Within this 400 million light-year volume of space are at least ~ 10 4 galaxies that 



possess a central bulge. Only a few percent of central galactic black holes at the current 
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epoch are AGNs or quasars. It is therefore plausible that most accretion onto massive and 
supermassive central galactic black holes systems at the current epoch is highly radiatively 
inefficient and dilute. 

Although the MRI is one of the more well-known free-energy gradient instabilities with 
astrophysical applications, in dilute astrophysical plasmas, such as radiatively inefficient 



accretion flows (RIAFs), large viscosities (Balbus 2004 Islam & Balbus 2005), and thermal 



conductivities (Balbus 2001) can also destabilize plasmas that possess free energy gradients. 



This unusual behavior arises from the fact that even a magnetic field with dynamically 
unimportant Lorentz forces can easily be strong enough that the ion Larmor radius is 
smaller than the collisional mean free path. Under these conditions, the viscous stress and 



thermal flux are dynamically important and directed along magnetic field lines (Braginskii 



1965). 



The explanations for the magnetothermal instability (MTI) and magnetoviscous 
instability (MVI) are conceptually simple. An equilibrium is perturbed so that there is 
a component of field line, hence of thermal conductivity (viscosity), along the thermal 
(angular velocity) gradient; heat (momentum) can then be transferred outwards, which 
further deforms magnetic field lines along the direction of the gradient - the process 
runs away. Similar behavior arises in the collisionless regime, in which large heat fluxes 
(momentum fluxes) along magnetic fields act to destabilize the plasma. This phenomenon 



is shown in collisionless and mildly collisional treatments of the MRI (Quataert et al. 2002 



Sharma et al. 2003). Furthermore, consistent collisionless MHD simulations (Sharma et al. 



2006, 2007) incorporate the fact that the arbitrary development of pressure anisotropy is 



limited by pitch-angle scattering due to fast gyrokinetic instabilities that violate magnetic 



adiabatic invariance for electrons and ions (Gary et al. 1994; Gary & Wang 1996 Gary 
et al.||1998l). 
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Numerical simulations have also begun to demonstrate how large anisotropic transport 
coefficients can destabilize a plasma that contains free energy gradients. More recent, 
comprehensive, local and global nonlinear simulations of the MTI in a standard MTI- 



unstable simulation box (McCourt et al. 2011) demonstrated steady states characterized by 



regions with relatively tangled magnetic field configurations, developed kinetic turbulence, 
driven by the free energy of an upwardly decreasing temperature profile. Furthermore, 



numerical nonlinear simulations of the collisionless MRI (Sharma et al. 2006 Sharma & 



Hammett 2006), which allows for efficient transport of angular momentum along field 
lines, also demonstrate the development of outwardly angular-momentum transporting 
turbulence, to a nearly steady state, in a local Keplerian rotational box (MRI-unstable, in 
that angular velocity decreases radially outwards). 

This paper expands on previous research by considering an incompressible fluid 
treatment of instability in a rotating plasma with large anisotropic ion viscous stresses as 
well as electron thermal conductivities. This is referred to as the magnetoviscous-thermal 
instability (MVTI). The organization of this paper is as follows. In §2 we demonstrate that 
the MVTI may play an important role in describing the dynamics of accretion flows about 
dim galactic nuclei, taking Sagittarius A* as an example. In §3 we lay down the equations 
of the system, describe our equilibrium state, and derive expressions for angular momentum 
and total energy balance in dilute, radiatively inefficient magnetized accretion flows for 
thin, largely rotationally supported disks. In §4 we derive the dispersion relation. In §5 
we justify and derive bulk fluxes for the MVTI, demonstrating that turbulence driven by 
this MHD instability can drive accretion onto radiatively inefficient flows. Finally, in §6 we 
provide a short discussion of the conclusions, its main astrophysical implications, as well as 
directions for further research; a companion paper will explore the MVTI in the collisionless 
limit, and the corresponding instability will be termed the collisionless MTI. 
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2. Preliminaries 

We briefly review the conditions of validity for the approximations used in this 
calculation, and note the forms of viscous stress and heat flux in a magnetized fluid. First, 
we note that a fluid approach for a low density plasma requires that the ion and electron 
mean free paths be small compared to global length scales. In our calculation, the ions and 
electrons are assumed to have the same temperature, so the ion-ion and electron-electron 
mean free paths are identical. For convenience, however, we will always refer to the ion 
mean free path Aj, 

T 2 

Xi = 1.5 x 10 13 — ±— cm, (1) 
rii In A 

where T 4 is the temperature in units of 10 4 K, n\ is the number density in cm -3 , and 
In A is the Coulomb logarithm. The second condition is that this mean free path be large 
compared to the ion Larmor radius Tl, 

r L = 9.5 x 10 7 £?^T 4 1/2 cm > ( 2 ) 

so that transport is directed along magnetic lines of force. Here B^q is the magnetic field 
strength in fiG. The electron Larmor radius is of course a factor of 43 smaller than the ion 
Larmor radius. 

The formal regime of validity is fx <C Aj <C R, where R is a typical global scale of the 
flow. In practice the lower bound on Aj is expected to be well-satisfied in the applications 
of interest (the protogalactic interstellar medium and under luminous black hole accretion). 
The upper bound on Aj, on the other hand, breaks down in the innermost regions of black 



hole accretion models for the Galactic Center (Quataert 2004). Our fluid approximation 
is expected to be well-satisfied throughout the bulk of this flow, from regions of a few 
times 10 4 Schwarzchild radii out to Sag. A*'s Bondi capture radius and the protogalactic 
interstellar medium. 
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3. Formulation of the Problem 

Dim mass-starved accretion flows onto black holes are radiatively inefficient. Under 
these circumstances, the turbulent energy flux should contain a substantial purely thermal 
contribution (from correlated fluctuations in the density and temperature) as well the usual 



contribution from the angular momentum transport (Balbus 2003). Because of its role 
in redistributing energy, this additional thermal flux could have a profound effect on the 
observational properties of the flow. It is the possibility that the MVTI could be the source 
of a profoundly enhanced thermal energy flux that motivates our study, beginning with an 
elucidation of its linear properties. 

We examine the rotational and convective stability of a disk under the effects of a 
magnetic field, under the effects of a large Braginskii viscosity ( Balbus|[2004 Islam & Balbus 



2005) as well as large thermal conductivity (Balbus 2001). Our coordinate system for the 



disk will be a standard cylindrical system: radius R, azimuth <ft, and axial variable z. As 
in 



Balbus (2001), we consider the stability of an equilibrium plasma at the midplane with 
only radial gradients in pressure and temperature, and nonradial equilibrium magnetic field 
B = B ^0 cos x4> + £ s i n x) ) where -Bo is the magnitude of the equilibrium magnetic field. 
We analyze unstable axisymmetric modes in the Boussinesq limit - incompressible flow and 
isobaric perturbations. We calculate the linear stability of an idealized plasma with a single 
temperature T with physical radially outwardly decreasing temperature and pressure that 
remains convectively stable, and demonstrate that its dispersion relation reduces to that of 
the MRI, MVI, and MTI in specific limiting cases. We also demonstrate that the quadratic 
heat fluxes and Reynolds stresses associated with the MVTI are of the right sense to drive 
accretion in radiatively inefficient astrophysical plasmas. 
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3.1. Constituent Equations 

The fluid equations consist of continuity, force balance, energy balance, and the 
induction equations: 

^ + V-( P V)=0, (3) 

fdV \ ( B-B\ BVB 

p + V • Vv = -V p + — — + — pV$ - V • cr, (4) 



V dt J \ An J An 

3 p dln^ = _ v _ (9b) _ <r:VV) (5) 

— =Vx(VxB), (6) 

where p is density; V is the velocity; $ is the gravitational potential; B is the magnetic 
field; cr is the viscous stress tensor; q is the anisotropic heat flux; and P is the pressure, 
b = ~Q/B is the unit vector along the magnetic field line. The viscous stress tensor and 



heat flux are given in Braginskii (1965): 



cr = -3pv (bh -^ij (h ■ VV • b - ■ , (7) 

q = -npb ■ \/9, (8) 

where v and k are the viscous and (electron) thermal diffusivity along the magnetic field 
line, 9 = p/p = k^T jm^ is the isothermal sound speed squared, and is the ion mass. 



3.2. Fluxes in Dilute, Magnetized Disks 

Here we derive expressions for the heat and angular momentum fluxes in geometrically 
thin, dilute, radiatively inefficient, magnetized disks. Since velocity fluctuations are possibly 
of the order of the disk sound speed, then the dominant contributions to turbulence that 
carries out angular momentum and (as we shall see) thermal energy to allow accretion is 
through quadratic correlations between fluctuating fluid quantities. We first find it useful 
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to transform to a frame corotating with the equilibrium flow V = RQ(R)(f), and define a 
flow velocity v = V — Vq with respect to the equilibrium. The equilibrium force balance 
equation is given by, 



„ /nl - 1 GM ~ GM A 

-R Q(R)R = Vpo - —^R - —^zz. 

Po R l R 6 



(9) 



If we employ a variable transformation V — > v = V — i?fi(i?)0, substitute in the forms of 
the viscous stress tensor (Eq. |7j) and heat flux (Eq. (8l), and take into account equilibrium 
force balance (Eq. [9]), then the dynamic MHD equations, Eqs. ^ - become, 



(i + n£|p+v.w = o, 



P 

B VB 



d d 
h — 



+ v ■ Vv - 2fiv xH Q'Ru R (p v = -V ( p + — J + 



£? 5 



•57T 



47T 



Vp - V • <T, 



' ^ fi^rVnpp 5/3 + -p\ • Vlnpp 5 3 



Of 



V • (gb) - cr : Vv, 



— + VL— ) B + v • VB = -B (V ■ v) + B • Vv + tt'RB R (f). 
at a 



(10) 

(11) 

(12) 
(13) 



We use a prescription for averaging out turbulent fluctuations as laid out in |Balbus & 



Hawley (1998). We assume a disk that is steady and spatially smooth over time and length 



scales much larger than the turbulence. The fluctuation in fluid velocities (u^, u z ), magnetic 
field B, and density p average to zero. 



(6114,) = (Su z ) 



0. 



(SB) t 



o. 



(8p) p = 0. 



(14) 
(15) 
(16) 



This average, () , is defined in the following manner: 



AT rR+AR/2 r 2w poo 



2tiT,ARAT J J R -ar/2 



Ap dR' dz d<fi dt, 



(17) 



J z=—oo 
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where, 



pdz. 



If 



This average is taken over a radial slice AR much smaller than the disk radius but larger 
than the largest scale of the turbulence and a time scale AT much longer than the turnover 
timescale of the turbulence. The radial velocity ur consists of a fluctuating part (Sur) = 
and a nonzero mean mass flow. The mean flow velocity is far smaller than the root mean 

( u b) p < \ (Su R 5u lj> ) f 



square turbulent velocity, i.e., 



The angular momentum density in the ^-direction, C z = pRV,p, can be found from 
Eqs. (|3| and Q. After some algebra, we have, 



dt 



Air 



P+ ^ ) ^ ,T " 



0. 



(19) 



If we substitute v for V and average according to the prescription given by Eq. (17), 
angular momentum balance is given by, 



The averaged quadratic angular momentum flux is, 

5Br6B, 



0. 



(20) 



(Tr^) = (oqSurSu^ - 



47T 



3pou [bo ■ V5v • b + Q'R5bR cos x] x <5&r cos x 



(21) 



The local total energy balance equation can be found in steps. First, if we dot Eq. (11) 



with v, and use Eq. (10), we have 



d d 

v • I — + fi— ) (pv) + v • Vp + Uidj (pUiUj) + 

Uidj ^—3pu 
P 



b . Vv • b V • v 

3 



bibj - 



+ Q'RpuRU^ 



(22) 



-v ■ Vpo — v ■ 



(V x B) x B 

47T 
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Second, if we dot Eq. (13) we get the equation for the magnetic energy density evolution, 

3\ B 2 



d_ 

dt 



8tt 



v V 



8tt 



— V • v H b • Vv • b+ 

47T 47T 



VL'R 



BrB 4 
Ait 



(23) 



Third, the MHD comoving internal energy equation, Eq. (12), is given by, 

° ' + ) P + V ■ f o^ v ) + pV • v = V • («p [b • V0] b) + 



2 Vd£ 



1. 



1, 



3pu b ■ Vv • b - -V ■ v b • Vv • b + Rh ■ VOZ^ - -V • v 



(24) 



We add Eqs (22), (23), and (24), we arrive at the total energy balance equation, 

B 2 ' 



d d 

h — 

dt d<t 



1 2 3 
2 PU + 2 P+ 8^ 



3pu 



b • Vv • b V • v 

3 



V- v 



1 



b(vb)--v 



1 2 5 
-pu 



1 



+ 



B x (v x B 



47T 



pv ■ — Vpo 

Po 



(25) 



dfl 
dlnR 



PUrU^ 



BrB 4 

47T 



3pu 



b ■ Vv • b + Rb ■ VQb^ - - V • v 

3 



bRk 



Applying the prescription for the average as given by Eq. (17) to Eq. (25), we arrive at the 
formula for the averaged total energy balance equation, 

dQ 



d{£) 1 d 1 dp 

+ - \P U R>, 



dt RdR 



po dR 



dlnR 



(Tr4>) ' 



(26) 



where, 



1 2 3 B 2 
2 PU + 2 P+ ^ 



(27) 



(F E r) = ^Po (Su R 59) p - Kp {5b R 5b R ^^ 



+ 



(28) 



Pqu ((b • V5v • b + Q'R5b R cos x) ^b R ) p . 
In Eq. (26) we have ignored the flux of gas kinetic energy since it appears as a cubic 



correlation in fluctuating quantities, and in Eqs. (21) and (26) we have ignored the Poynting 



flux B x (v x B) / (47r) since it is subdominant to other terms in the angular momentum 
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and energy flux. We have ignored terms of the form V ■ v in the viscous stress, which 



appear in Eqs. (21) and (26), since we are in the Boussinesq limit. In the absence of a 
radiative channel by which the power generated through accretion, — (dQ/dlnR) (Tr^), is 
dissipated or radiated, the energy carried away by the turbulence {Fer) must be positive. 



3.3. Dispersion Relation 

Assume small perturbations 5a about equilibrium fluid quantities ao; thus v = 8v, 
B = B + <5B, 6 = 9 + 59, and p = po + 5p. We assume nonaxisymmetric perturbations of 
the form: 

5a oc exp (ikjiR + ik z z + "ft) . (29) 

In equilibrium the viscous stress and heat flux are zero. One can demonstrate that the 
perturbed viscous stress tensor and heat flux are given by the following, where <5B = 5~Q/ Bq\ 



(30) 



5a = -3p u (bb - ^ij (i (k ■ b) (6v ■ b) + Q!R5B R cos X ) , 

dq = - K p dB R —Q^ Kp (ik-b)— . (31) 

The perturbed induction equation, from Eq. ^ reduces to the following: 

j5B R = ik z sin x5v R 

"fSBfi = ik z sin x$ v 4> + &'R5Br (32) 
-f5B z = ik z B sin x5v z , 



From Eq. (32) the perturbed viscous stress tensor has the following nonzero components: 

5a H = 5o zz = ^cos 2 X - <5°bb, (33) 

5cr RR = --5a hh , (34) 
5a Z(P = 5a^ z = sin \ cos x^bb, (35) 
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where 5o"bb is the anisotropic pressure, or equivalently <5<7bb/2 is the component of viscous 
stress along the magnetic field: 



The perturbed energy balance equation, from Eq. ^ is given by the following: 

<91np Po 5/3 



(36) 



5p 5 Sp 
Po 3p 



Sv R - 



dR 



k ik z smx oB R 



OR 



1 2 • 2 
- K sin x 



66 

Tn 



x Sa : V (m$\ . 



(37) 



First, the form of the perturbed stress tensor as given in Eqs. (33), (34), and (35) implies 
that Scr : V (^RflcftJ = 0. Second, Eq. (37), with isobaric perturbations, implies the 
following perturbed density: 



Sp 3 Sv R d\np p 5 ^ 3 /dR — ^ink z smxSB R d\nT /dR 
Po 5 7 + \nk\ sin 2 % 



(38) 



Therefore, the perturbed form of the force balance equation, Eq. Q 



dSv 



20,% x <5v + VL'R4> 



-V Sp 



Po 

5p e u dlnp ° R 



h-SB 

47T 



b • V5B 



4vrp c 



p X V ■ Sa+ 



(39) 



Po dR 
In component form this reduces to, 

' 5p B cos x^Bcf, + B sin xSB z 



jSvr - 2£lSv 4 



ik R 



Po 



4vrp 



+ — oB R - 

4vrp 



ik R p da RR -ik z p Sa zR -\ ^0 

Po OK 



j5v (f , + (2Q + Q'R) 5v R = ikz ^° SmX §B (j) - ik R p l Sa <t>R - ik z p 1 Sa 4 



jSv z = - ik z 



Arcpo 

5p B cos x^Btf, + Bq sin xSB z 
Po 47rp 



ik z B sinx 
47rp 



5 Br- 



im 

(41) 
(42) 



ik R p 1 Sa zR - ik z p L Sa 



-1 



Eq. (42) results in the following expression for the total pressure, where v\ = Bq/ (4irp): 



Sp 
Po 



+ v\ cos x^Bif, =wy (3 sin 2 x - l) ( £-B</> cos x - 'tt&B r sin x ) - -t£i& v R, (43) 
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From the induction equation, Eq. 32 , the force balance equations, Eqs. (40, 41, 43) can be 



rewritten in terms of 5B R and 5Bs'- 



k 



1 + j| J 1 2 8B R - 2Q (jSB^ - n'RSBn) = 3vjk R k z sin 3 X [ 5B<j> cos X -f8B R ^ X \- 

(44) 



{k 2 R + k 2 z ) v\ sin 2 X 5B R +-9 



dR 



3 n fd\n Po \ ^d\np p- 5/3 /dR+ \^k 2 z sin 2 x<9 In T /dR 



2,_jr,2 • 2 



7 + \^k 2 z sin 2 X 
k R 



5B 



R, 



7 dBft, + 2Vt^5B R = —3i"yk z sin X cos x I ^-B^ cos X — —5B R sin x I — k z v A sin X 5B^. (45) 



k? 



This results in the following dispersion relation: 



k 2 2 dn 2 

— 7 H 

&2 (ilni? 



+ 3vk R ^i sin x — ?>vk R k z VL' R sin x cos x + & sin x — 



3„ /<91np \ 7^1np Po 5/3 /^+ 1«^ 2 sin 2 x<91nT /&R ' 



:0n 



dR 



7 + \^k\ sin 2 x 



X 



(46) 



( 7 2 + k 2 z v 2 A sin 2 x + 3^fc 2 7 sin 2 X cos 2 x) + 7 2 (4fi 2 - §v 2 k 2 R k 2 z sin 6 x cos 2 x) = 0. 
The Prandtl number for our regime of interest, where electrons dominate the thermal 



conductivity, is given in Braginskii (1965): 



Pr = v/ n 



0.96 /2m. 



1/2 



1/101. 



(47) 



3.2 \ rrii 

After making the following normalizations, where ctp, as, and «t are the normalized inverse 
scale heights of pressure, entropy, and temperature, respectively: 

k = \<.va/Q 

7 = j/Q 
v = vVl/v 2 a 



k = kVL/v 2 a 

d In p 



(48) 



ap = —H 
o:t = —H 

a s = -H 



dR 
d In T 

OR 

din p p~ 5/3 _ 5 
OR ~~ 3 



OCT - ^OSP, 



- 15 - 



i /2 j y 

where H = 9 /ft is the disk scale height. One can show that Eq. (46) reduces to the 
following: 



-f + 2 



; - • '• • /, - Mil 2 \ - 



dlnR 



+ 3z>/c R 7sin \ ~ 3vkpk z 



d\nR 



3 cts7 + f a^Pr ^fc 2 sin" 1 x 
-ap- 



7 + |Pr 1 vk 2 sin 2 x 
7 2 (4 - 90 2 k 2 R k 2 z sin 6 x cos 2 x 



x I 7 + A; z sin x + 3^k 2 7 sin x cos X ) + 



0. 



(49) 



Physically speaking, in this fluid analysis we are always in the regime in which k > u, hence 
where Pr > 1. The result is that for sufficiently large transport coefficient vVl/v 2 A > 1 we 
then have a density perturbation given by: 

6p 5vr ( <91nT" 



A) 



7 



OR 



and the following substitution occurs in Eq. (49): 



3 a si + farPr 1 vk 2 sm 2 x 
-a P —7- > a P a T . 



7 + |Pr Pk^sin 2 x 



Note the limits of Eq. (49): 1) in the limit that the transport coefficients go to zero and 
there are no equilibrium gradients, we reproduce the MRI; 2) in the limit of no equilibrium 
gradients and finite viscosity, we reproduce the MVI; 3) in the limit of no viscosity but 



finite thermal conductivity, we reproduce the MTI dispersion relation (Balbus 2001); and 
4) in the limit of no thermal conductivity or viscosity, but finite equilibrium gradients, we 
reproduce the dispersion relation for convectively unstable modes in a rotating magnetized 



plasma (Balbus 1995). To put these parameters in perspective, in physical terms the 
condition under which the magnetoviscous instability operates is one in which uQ/v A > 1. 



The viscosity and Alfven velocity in dimensional units is given by the following: 



[/ = 1.4x 10 



l!) -M 



5/2 



cm 2 s , 



Va 



Hi In A 

2.2 x 10 5 B^gtii 1/2 cm s -1 . 



(50) 
(51) 
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One must thus have an angular velocity of the flow such that: 

n > 3.5 x 1(T 9 B\ G (lnA)T 4 ~ 5/2 s _1 , (52) 

for the viscosity to be dynamically important. Furthermore, the thermal diffusivity k ^> v, 
so that the range of applicability of thermal effects is significantly larger than that of 
viscous transport. 



4. Results 

In this section we calculate the quadratic azimuthal stress (Tr^) and radial heat flux 
qu associated with these instabilities. Since the MVTI operates in a regime of high thermal 
diffusivity and possibly high viscous diffusivity we need to consider the local energetics and 



angular momentum transport modified by large viscous and thermal fluxes. From Eqs. (21) 



and (28), we estimate the quadratic angular momentum flux and heat flux carried away by 
a mode of wavenumber k as the following, 

(T R<t> ) = Re (po6u R 6ul - p v 2 A SB R SBl- 
"iip^v [(k ■ bo) (5vk • bo) + Q'R5B R cos x] $B*r cos X) 
(F E r) = Re C^po6u R 68£ - Kp ^\SB R \ 2 + 

ip v [(b • Sv k ) (b • k) + Q'R5B R cosx] 5B* R 

8a refers to fluctuating quantity a for an MVTI eigenmode of wavenumber k. 



(53) 



(54) 



4.1. Growth Rate and Stability Characteristics 



Here we consider two equilibrium configurations that illustrate the magnetoviscous- 
thermal dispersion relation. We only consider a physical pressure profile, hence a R = 10 
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(i.e., outwardly decreasing pressure). We also consider a flow that is convectively stable, 
hence with HdlnPp~ 5 ^ 3 /OR > or q§ < 0. This implies the following: 



< a? < -Op, 
5 

which implies that ay < 4. We focus our results on only vertical wavenumbers, 



(55) 



kpi = 0. Then Eq. (49) reduces to the following quintic polynomial, where we denote 



x = (k • va) /SI = k z sin x- 



9 dlnQ 9 



7 H — Pr 1 z>x z 
5 



«s7 + -a^Pr x i>x 2 



x 



(56) 



(7 2 + x 2 + 3z>7x 2 cos 2 x) + 47 2 ^7 + jjPr^z/x 2 ^ = 0, 

In Fig. §, we examine the instance where both magnetic tension and thermal conductivity 
are dynamically significant; dispersion relations match the salient characteristic of the 



magnetoviscous instability (Balbus 2004 Islam & Balbus 2005) - saturation of the growth 
rate at small wavenumbers, described by the condition vk 2 ~ Q. For the case «y = 
(no equilibrium temperature gradient) we approach the magnetoviscous instability - the 
mode is reproduced only for ap = and «t = 0. In Fig. g), we examine the instance 
where only the thermal conductivity is dynamically important; the dispersion relations 
are similar to the magnetorotational instability. From an examination of the dispersion 



relation Eq. (46) or the normalized dispersion relation Eq. (56) that for physical ax > the 
growth rate and wavenumber extent of the instability are increased - this is shown also in 
Fig. Q and Fig. One expects that equilibrium scale heights of temperature, pressure, 
and entropy are of order the radius R. Thus, in order to have significant magnetoviscous 
and magnetothermal effects one requires that a^, ap, and as be at least of order unity. 
Hence, only thick disks H ~ R are expected to be susceptible to these classes of instability. 
However, note that we have calculated quadratic modal fluxes of ostensibly thick disks using 
thin-disk expressions for angular momentum (Eq. [53]) and heat (Eq. [54]). This may not 
pose an intractable problem: first, even for relatively slender disks, ones in which H ^ R, 
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01234567 

Fig. 1. — Plot of the real portion of the growth rate for various olt for a Keplerian-like 
rotation profile, VL oc i? _3//2 , significant viscous diffusion coefficient vVL/v 2 A = 10 2 and Prandtl 
number Pr = 1/101 (see Eq. |47j ) - this corresponds to the case where both anisotropic 
viscous forces dominate over magnetic tension. 

we can expect that the angular momentum and heat fluxes due to the nonlinear MVTI will 
still be of the same sign (outwards) and magnitude as their respective thin-disk expressions; 
and second, even for thin RIAFs, the nonlinear MVTI or its collisionless analogue appear to 
provide the only plausible MHD mechanism for the transport of accretion energy outwards. 
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k-v A /n 

Fig. 2. — Plot of the dispersion relation for various ay for a Keplerian-like rotation profile, 
oc R~ 3 / 2 , a small viscous diffusion coefficient vVt/v\ = 1 and Prandtl number Pr = 1/101 
(see Eq. [47]). 

4.2. Heat Flux and Azimuthal Stresses 



The following are the calculated perturbed velocities, magnetic fields, density, and 
temperature in terms of the radial displacement £r = ■y~ 1 5vR, that employing the induction 
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equations (Eq. (32] ) and Eqs. (38) and (45). We examine nonradial modes, k R /k z = 0: 



5v R = 7 ^£r 
5v a = - 



if 



7 2 + x 2 + 3z>x 2 7 cos 2 x 
5B R = i(k-B)£ R 

2 7 



+ 



dlnft 
dlnR 



B 



5B, 
5p 



-i (k- B) 



B 



(57) 



7 2 + x 2 + 3i>x 2 7 cos 2 x 
3 7a 5 + |Pr _1 z>x 2 a T - , 



7 + §Pr 1 Pa; 2 



3 7"5 + |Pr l vx 2 a T 
-OoSp/po = - x — ; f- 

5 7 + fPr Vr 



with the form of the perturbed anisotropic pressure from Eq. (36), the turbulent heat flux 



from Eq. (54) yields the following: 



Feb. = PoRe ( -S95v* R + p wy cos x$v R 5B, 

3 A ia s + pPr~ l x 2 a T i ,-, , 

7^7 . , 2-p -i 2 ^ # 

2 7 + |i/Pr or 



PoPr Vfc z Im (sin x<5#51?^.) 



(5? 



And the azimuthal stress from Eq. (53) yields the following: 



SBrSBI _ 
Re PqSvrSv^ — 3p z/cos xi^B* R 5B q 



A.TX 



7 2- 



4f 



(59) 



p n 2 |e 



dlni? 7 2 + x 2 + 3z>x 2 7 cos 2 x. 
In all the quasilinear plots we keep ap — 10, choose > for which the flow is convectively 
stable, and fix the Prandtl number Pr = 1/101. In Figs. ^ and Q are plots of the 
normalized thermal flux q R and azimuthal stress T R( p for magnetoviscous-like dispersion 
modes [yVl/v 2 A = 10 2 , see Fig. [I]). In Figs. ^ and ^ are plots of normalized q R and 
T R( p for magnetorotational-like modes {vVt/v 2 A = 1, see Fig. |2j). From Figs. ^ and ([5]), 
one observes that the quadratic heat flux in the limit of large {yVl/v 2 A S> 1) and moderate 
(uQ/v\ ~ 1) viscosities results in dynamically significant outwards heat fluxes, with 



-21- 




01234567 

k-v A /n 

Fig. 3. — The normalized radial flux of thermal energy with parameters as described in 



Fig. Q), where uVt/v\ = 10 2 . 

qn ^ Tjitj, (kBOo/rrii) 1 ^ 2 . This implies that the nonlinear MVTI can play an important role 
in transporting out the energy, generated via the azimuthal stress, in nonradiative accreting 
flows. Second, from Figs. Q andm ([6]), one notes that the flux of angular momentum for 
these modes can be either inwards or outwards. One also observes for a range of small 
wavenumbers that thermal gradients may cause a negative azimuthal stress, unlike the 
MRI. 

An important point that must be made is that with a large thermal conductivity, 
angular momentum momentum may be transported inwards or outwards for large decreasing 
gradients of temperature. This effect can be seen even for the purely magnetothermal 
instability (MTI) in a rigidly rotating plasma. We can re-derive the dispersion relation and 
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Fig. 4. — Normalized azimuthal momentum radial flux with parameters as described in 
Fig. 0, where vVt/v\ = 10 2 . 

angular momentum flux for the MTI by setting Pr _1 — k/u, while letting v — > 0. Then the 



normalized dispersion relation, Eq. (56), and azimuthal stress, Eq. (59): 

2 



r + 2 



dlnfi 
dlaR 



7 + -kx 



:Otp 



a Sl + -Ot T KX 

O 



(f + x 2 ) + 



Atf ( 7 H — kx 2 
5 



^ = 7 



x 2 (2-dlnft/dln.R) - 7 2 (2 + dlnft/d]n.R) 



7 2 + x 2 



R\ 



(60) 



(61) 



Fig. ([7]) demonstrates that in a rigidly rotating plasma = 0, the magnetothermal 
instability can drive a positive azimuthal stress depending on wavenumber. As the system 
approaches marginal convective stability as — >■ from isothermality (ay = 0), the range 
of wavenumbers for which the stress is outwards decreases. One must note, however, that 
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Fig. 5. — The normalized radial flux of thermal energy for a magnetorotational-like mode, 
with parameters as described in Fig. (|2|. 
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Fig. 6. — Normalized (Tr^) for a magnetorotational-like mode, with parameters as described 
in Figs. @. For magnetoviscous-like modes, it is apparent that there is a range of low 
wavenumbers for which (Tr^) is inwards. 
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in the absence of rotational shear no energy can be extracted from the flow. Second, the 
ambiguity of angular momentum transport for the MVTI in this case is not necessarily 
evidence of unusual physics: the MVTI acts as a mechanism to transport thermal energy 
outwards, which is largely (and, in the case of rigid rotation, completely) independent of 
the manner in which it transports angular momentum. A first step into a numerical study 
to characterize MVTI stability could begin with a high spatial grid resolution simulation, at 
especially large Reynolds numbers, about a near-MVTI stable rotating system, in order to 
discover unambiguous (not numerical diflusivity limited) levels of turbulence and angular 



momentum and heat flux. This would be analogous to Lesur & Longaretti (2005), who 
have demonstrated the unambiguous presence of minuscule turbulence in high Reynolds 
number (Re ~ 10 5 ) hydrodynamic numerical simulations of differentially rotating flows near 
Rayleigh stability. 
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k-v A /n 

Fig. 7. — Normalized flux (Tr^) / (pf2 2 |£r| 2 ) for a rigid rotation profile (Q'R = 0) for the 
magnetothermal instability (see Eq. [60j). There is a much larger range of negative (Tr^) 
(inwards flux of angular momentum) than shown in Fig. ([6]). 

5. Conclusions and Further Research 

In this paper we attempted to characterize the stability of a mildly dilute rotating 
accretion flow to axisymmetric magnetoviscous-thermal modes, by incorporating both an 
anisotropic magnetized viscous stress tensor and anisotropic electron thermal conductivities. 
We have demonstrated several important properties of the MVTI. First, we find that thick 
disks with sufficiently large thermal diffusion coefficients are susceptible to these classes 
of temperature gradient instabilities. Second, physical outwardly decreasing temperature 
gradients increase the range of wavenumbers for which the magnetorotational-like instability 
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(where the viscous diffusion coefficient is small vfl < v\) and magnetoviscous-like instability 
(where the viscous diffusion coefficient is large) are unstable, as well as increasing the 
growth rate of unstable wavenumbers. Third, quadratic fluxes of angular momentum are 
substantially modified by both a large outward viscous transport in dilute plasmas. Fourth, 
in rotating systems with equilibrium outwardly decreasing temperature, the linear MVTI 
can transport angular momentum either inwards or outwards. 

Companion papers will consider the following. First, we will analyze a more general 
case of the collisionless and mildly collisional MTI, demonstrating that there is an effective 
viscous and thermal diffusivity of order Oo/Vt along the magnetic field lines. Second, we 
demonstrate the dynamics of the collisionless and mildly collisional MTI with electron 
pressure dynamics, in order to characterize the level of electron (and ion) heating due 
to nonlinear collisionless MTI turbulence. And third, we consider more unusual physical 
phenomena (finite equilibrium acceleration of ions and electrons away from the disk 
midplane, or collisionless reconnecting modes in a Keplerian profile) in collisionless disks. 
A necessary next step would be to consider more carefully an appropriate local equilibrium 



profile (see, e.g., Bisnovatyi-Kogan & Blinnikov (1972); Bisnovatyi-Kogan & Seidov (1985)). 



Our analysis of the MVTI may allow for a clearer physical resolution of analytic models, 



such as proposed advection dominated accretion flows (Narayan & Yi 1994), convectively 



dominated accretion flows (Quataert & Gruzinov 2000), or adiabatic inflow-outflow 



solutions (Blandford & Begelman 1999), of largely rotationally supported RIAFs. Energy 



transport for Rayleigh-stable flows is generally outwards where (dlnT/dlnR) > 0, and 
the maximal rate of quadratic energy transport Fer goes as H/R(d\nT/dlnR). This 
implies that the level of MHD turbulence and the level of energy transport in RIAFs may 
be some function of this parameter. In analytic models of accretion in dilute rotationally 
supported flows, a more appropriate estimate of the turbulent energy transport flux may be 



-28- 



H/R(d\nT/d\nR)p9 1 ' 2 ] this is in contrast to an energy flux estimate of pO 1 ^ 2 , applicable 



to largely pressure supported or spherically symmetric flows (Tanaka & Menou 2006 



Johnson & Quataert 2007). For geometrically thin RIAFs, due to quadratic estimates 



dependent on the parameter Hd\nT/dR, our results suggest that accretion will be far 
below the gravitational capture rate of matter from the black hole's ambient medium. 

Although much work remains to be done in understanding the dynamics of hot dilute 



accretion flows, this paper and its companions (Balbus 2001, 2004; Islam & Balbus 2005) 



have demonstrated a class of instabilities that may roughly describe the type of turbulence 
expected within RIAFs. Unfortunately, MHD simulations of a local radial slice of a 
Keplerian disk are ill-suited to studies of phenomena in which both angular momentum and 
energy transport must be characterized, as the total energy within the radial slice increases 
at a rate — (dQ /dhiR) (Tr^) . Global collisionless MHD simulations, employing fluid codes 



that naturally conserve total mechanical energy such as ATHENA (Stone et al. 2008), show 



the most promise in beginning our understanding of the nonlinear stages of these newly 
explored free-energy channels, and possibly in demonstrating accretion via nonlinear MVTI 
turbulence. 
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